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ABSTRACT 

We take advantage of more than 10 years of monitoring of the eclipsing HMXB systems EMC X-4, Cen X-3, 4U 1700-377, 4U 
1538-522, SMC X-1, IGR J18027-2016, Vela X-1,IGR 117252-3616, XTE 11855-026, and OAO 1657-415 with the ASM on-board 
RXTE and ISGRI on-board INTEGRAL to update their ephemeris. These results are used to refine previous measurements of the 
orbital period decay of all sources (where available) and provide the first accurate values of the apsidal advance in Vela X-1 and 
4U 1538-522. Updated values for the masses of the neutron stars hosted in the ten HMXBs are also provided, as well as the long-term 
lightcurves folded on the sources best determined orbital parameters. These lightcurves reveal complex eclipse ingresses and egresses, 
that are understood mostly as being due to the presence of accretion wakes. The results reported in this paper constitute a database to 
be used for population and evolutionary studies of HMXBs, as well as theoretical modelling of long-term accretion in wind-fed X-ray 
binaries. 
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1. Introduction 


High-mass X-ray binaries (HMXB) are among the brightest X- 
ray sources in our Galaxy, and were discovered for the first 
time with the Uhuru satellite ( [Giacconi et al.|[r971 1 and inde¬ 
pendently with balloon-born instruments (see e.g., Lewin et al. 
|1971| l. These sources typically host a compact object (often a 
neutron star, NS) accreting material lost by a massive compan¬ 
ion. Depending on the nature of the latter star, HMXBs are di¬ 
vided into Be or Supergiant X-ray binaries (SgXBs). In the first 
case, a conspicuous X-ray emission is generated due to the ac¬ 
cretion of the matter from the dense decretion disk around the 
companion star, that leads to transient events (outbursts) of dif¬ 
ferent types (see e.g., Reig||20lT for a review). In the SgXBs, 
the compact object can either accrete through the fast and dense 
wind produced by the companion, or in a few cases through the 
so-called Roche-Iobe overflow (a number of systems showed ev¬ 
idence for both mechanisms contributing to the accretion onto 
the compact object and the overall X-ray emission; see, e.g., 
|Chaty|20TT| for a recent review). 

In Roche-lobe accreting systems, the inflowing material 
from the companion star is endowed with a high specific an¬ 
gular momentum, and thus accretion is generally mediated by 
a disk surrounding the compact object. As this is a very effi¬ 
cient way to accrete matter, disk-fed systems achieve usually a 
relatively high X-ray luminosity (~ 10^^ erg s“^). In wind ac¬ 
creting systems, the compact object is instead deeply embedded 
in the material lost by the companion star, which do not posses 
sufficient angular momentum to form an accretion disk and ac¬ 
cretes quasi-spherically onto the NS (see, e.g. Bozzo et al.|2008[ 


Shakura et al. 2012 for recent reviews). This process lead to 


a strongly variable X-ray emission, reaching values significantly 


lower than those attained by dis k-accreting system (ranging from 
lO^'* to 10^® erg s“*; see, e.g., 

|T99^|Bildsten et al.|19971 i. 

Since their discovery, HMXBs have been intensively mon- 


Joss & Rappaport||l984 Nagase 


itored in order to test wind accretion models onto magnetized 

NSs (see e.g., |Lamers et al. |1976[ 

Bhattacharya & van den 

|Heuvel||1991[ |Lutovinov et al.||2013 

and references therein). 


measure the long term spin-variation of NSs as function of the 
mass accretion rate (probing the so-called “accretion-torques”, 
see e.g., Ghosh & Lamb|1979[ Lovelace et al.|1995) l, and study 
the evolution of the orbital parameters of these system through, 
e.g. measurements of the orbital period decay (see|Bildsten et al.| 


1997 and references therein). In the majority of HMXBs, pul¬ 


sations in the X-ray emission firmly established the presence 
of NSs as accreting compact objects (see e.g., |Lutovinov et al.| 
2005 Lutovinov & Tsygankov 2009 1 for a recent review and 
HMXBs statistic), and the detection of cyclotron absorption 


netic field strength (~ 1 - 5 x lO'^ G; see, e.g.. 

Cobum et al. 

2002 Filippova et al. 20051 |Caballero & Wilms| 

2012|l. In the 


ence of black holes companions cannot still be firmly ruled out 
(see, e.g., the case of 4U 1700-377, |Rubin et al.|1996[ Clark et al. 
| 2002 | ). 


Among more than hundred HMXBs, only in a few sources 
the inclination of the system is high enough that the compact star 
is periodically occulted along our line of sight by the compan¬ 
ion, giving rise to X-ray eclipses. For these sources, it is possible 
to infer a number of orbital parameters (e.g., orbital period) from 
the measured duration of the eclipse. The energy-dependent pro¬ 
file of the X-ray lightcurve during the eclipse ingress and egress 
also reveal details of the OB stellar wind structure (e.g., |White| 
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Table 1. Orbital parameters of the 10 eclipsing HMXB systems in ascending order of the orbital period. Where relevant, we indicated 
in brackets the uncertainties on the last digits of each reported value. 


Source 

Type 

Orbital Epoch'" 

PoYb 

P. 

P orb / P orb 

flx sin i 

e 

OJ 



(MID) 

(days) 

(sec) 

10-" yr-' 

(It-sec) 


(deg) 

LMC X-4“ 

08III 

51110.86579(10) 

1.40839776(26) 

13.5 

-0.98(7) 

26.343(16) 

0.006(2) 

- 

Cen X-3'’ 

0(6-7)11-111 

50506.788423(7) 

2.087113936(7) 

4.8 

-1.799(2) 

39.6612(9) 

< 0.0016 

- 

4U 1700-3//" 

06.51af-i 

48900.373(15) 

3.411581(27) 

- 

-3.3(6) 

48-82 

-- 

49(11) 

4U 1538-522'' 

B0.21a 

52851.33(1) 

3.728382(11) 

526.8 

- 

53.1(1.5) 

0.18(1) 

40(12) 

SMC X-L 

BOsg 

50324.691861(8) 

3.89220909(4) 

0.71 

-3.402(7) 

53.4876(9) 

0.00089(6) 

166(12) 

SAX 11802.7-2017' 

Bllb 

52168.26(4) 

4.5696(9) 

139.6 

- 

68(1) 

- 

- 

XTb, 11855-026*' 

BOlaep 

51495.25(2) 

6.0724(9) 

360.7 

- 

80.5(1.4) 

0.04(2) 

226(15) 

Vela X-l“ 

B'0'.51ae 

48895.2186(12) 

8.964368(40) 

283.2 

- 

113.89(13) 

0.0898(12) 

152.59(92) 

HXO 1722-363" 

BO-Blla 

53761.68(4) 

9.7403(4) 

413.9 

- 

101(1) 

<0.19 

- 

"OAO 1657-'415' — 

B0-6sg 

50689.116(50) 

10.44749(55) 

37.3 

-3.40(15) 

106.10(2) 

0.1033(6) 

87.6(1.3) 


Notes.l^ evine et al. 1 2000 1 ; Bildsten et al.|l 1997t;^bildst en et al.|( 1997t; Raichur & Paul j2010|;l|l^bin et al. 1 1996 1 ; Hammerschlag-Hensberg e 


et al. 


et al.. . I_I, I_^ I_I i-t^^=^=i-i=ii‘-te=_i=^=^^=i-t= 

^or )et & Mukai 1 2002 i; l^a.k.' a. IGR J17252-3616) Zurita Heras et al. (2006 1 ; Thompson et al. 


(lioOSk^Mukheriee et al.|^2006^ ; |C^ark|^2b00^;lfinW et al.ll20lb1i;|Raichur & PaulM2oio>;Ra.k a IGR J18027-2016) 


i 20 ll^ ; |Jenkeend:| ( 12011 


.I3rhe eccentricity for this source has been reported by 


Hill et al. 


|2TO7 Manousakis & Walter 1 2011 1 : 


Hammerschlag-Hensberge et al. 1 2003 i to be 0.22(4) and 


12005 i: 


Baykal 


later questioned by|Clark| 2000 1 and this work.l^lVIid-eclipse time, equivalent to time when mean longitude I equals 7t/2 for a circular orbit. 


|et al.|1995| ). In those systems in which orbital ephemerides can 
be accurately measured over several decades, it is expected that 
changes in the orbital period can also be revealed, providing fur¬ 
ther insight into tidal interaction and mass transfer mechanisms 
that regulated the accretion onto the compact object. In eccen¬ 
tric systems, accurate ephemerides can also be used to infer the 
angle of the periastron and the apsidal advance (see e.g., |van der| 
|Klis & Bonnet-Bidaud|1984| and references therein). 

In this work we take advantage of the available long-term 
monitoring observations of ten eclipsing HMXB^ carried out 
with the hard X-ray imager IBIS/ISGRI (17 keV - 200 keV) 
on-board INTEGRAL and the All Sky Monitor (ASM) on¬ 
board RXTE (2-12 keV), to obtain the most accurate ephemeris 
available so far for these sources ( Liu et ah] 2006| l. We sig¬ 
nificantly improved their previously measured orbital periods, 
eclipse durations, and masses of the accreting NSs hosted in 
them. Measurements of the orbital period decay and apsidal mo¬ 
tion is presented for the first time for several of these systems, 
and refined measurements are provided for all the others with re¬ 
spect to those already available in the literature. We discuss the 
energy-dependent profile of the long-term lightcurve during the 
eclipse ingress and egress for all these systems in terms of X-ray 
absorption expected in wind and/or disk accreting systems. 

1.1. The 10 eclipsing HMXBs 

In Table we list the ten sources selected for the present work 
and provide a summary of their main characteristics, as mea¬ 
sured in the most recently available literature: the orbital epoch, 
the orbital period Horb, the spin period P^, the orbital period 
change Porb/Porb, the projected semi-major axis sin i, the ec¬ 
centricity e, and the periastron angle o). All these systems show 
evidence for alternate spin-up and spin-down phases. In some 
cases the average spin evolution led to an effective long term in¬ 
crease or decrease of the NS spin rate (see e.g., [Bildsten et al.| 
1997| Inam et al.|2010 for more details), and thus the spin peri¬ 
ods reported in Table should only be considered as indicative. 


The spin period of the NS (possibly) hosted in 4U1700-377 is 
currently unknown. 

For the two HMXBs EXO 1722-363 and SAX J1802.7- 
2017, we indicated in Table [T] also the associated INTEGRAL 
source name IGR J17252-3616 and IGR J18027-2016, respec- 


' Be X-ray binaries are not expected to show long detectable X-ray 
eclipses, as the radial extent of the companion is less than those of 
supergiant stars in HMXBs and the orbital periods are usually much 
longer and characterized by high eccentricities i Reig 20111. 


lively ( Zurita Heras et aL]|2006[[Revnivtsev et al. 20041. These 
two sources are classified as highly absorbed HMXBs (see e.g., 
Walter et al. 2006| l, and indeed these are the only sources in 
which no eclipse is evident from their soft X-ray lightcurves 
(<12 keV, see Sec.|^. 

The 10 selected sources comprise both disk and wind accret¬ 
ing systems. A different shape of the eclipse ingress and egress 
is expected in these two cases as a function of the energy (see 
Sec.|^. The relatively large X-ray luminosity of LMC X^, Cen 
X-3, and SMC X-1 makes these sources the prime candidates 
for being disk-fed systems whereas 4U 1700-377, 4U 1538- 
522, SAX J1802.7-2017, XTE J1855-026, Vela X-1, and EXO 
1722-363 are all thought to be wind-fed systems. OAO 1657- 
415 is unique among the known HMXBs as it is believed to be 
a wind-fed system for most of the time, and undergo only spo¬ 
radically episodes of accretion from a temporary disk. Since the 
binary is too wide for Roche lobe overflow to occur, this may 
provide the first clear evidence that winds in HMXBs possess 
sufficient angular momentum to form accretion disks (see e.g., 
Bildsten et al.|1997| Chakrabarty et al.|1993| l. Eurther additional 
information on the ten sources can be found in [Liu et ari(|2006|l. 


2. Observations and data 

All the eclipsing HMXBs reported in Table [T] have been contin¬ 
uously monitored in the X-ray domain by the ASM since the be¬ 
ginning of 1996 and by the IBIS/ISGRI since the early 2003. We 
used public available IBIS/ISGRI lightcurves retrieved from the 
on-line tool Hi^-Energy Astrophysics Virtually ENlightened 
Sky {HEAVEN^^ HEAVENS data reduction was performed us¬ 
ing the standard Offline Science Analysis (OSA) version 9.0 dis¬ 
tributed by the INTEGRAL Science Data Center ( jCourvoisierj 
et al. 20031. Eor each source we downloaded the ISGRI high 
energy lightcurve binned over each pointing (science window) 


http://www.isdc.unige.ch/heavens/ 
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of roughly 2 ks, in the 17-40 keV and 40-150 keV band. The 
RXTE/ASM lightcurves have been retrieved from the NASA 
HEARS AC FTP servet]^ binned dwell-by-dwell (90 s bins) in 
the 1.5-3 keV, 3-5 keV, and 5-12 keV energy bands. In Table 
we report the observation time interval and the total effective ex¬ 
posure time of the ASM and ISGRI data for each of the selected 
source. All the photons arrival times of the INTEGRAL/lSGRl 
and RXTE/ASM lightcurves were corrected to the barycenter 
reference time of the solar system. We used the OSA task 
barycent for the ISGRI data, and the procedure described in 
the RXTE cookbool^ for the ASM data. The barycentric cor¬ 
rection is usually considered mainly to perform highly precision 
timing analysis on a short observational time interval. However, 
in the present case, it was required in order to ensure uniformity 
over a long-term set of data spanning more than a decade. 

For each of the eclipse found in the ISGRI and ASM 
lightcurves we measured the ingress, egress and the mid-eclipse 
time as described in Sec. 13.11 


Table 2. Log of INTEGRAL/lSGRl and RXTE/ASM observa¬ 
tions of the 10 selected eclipsing HMXBs. 


Source 

ISGRI 

(MJD-50000) 

Exp. 

(ks) 

ASM 

(MID-50000) 

Exp. 

(ks) 

LMC X-4 

2641.40-5005.48 

659 

87.30-5748.11 

5627 

Cen X-3 

2650.54-5155.90 

2033 

87.29-5747.29 

5297 

4U 1700-377 

2668.25-5256.35 

2863 

88.11-5749.81 

4709 

4U 1538-522 

2650.74-5256.37 

2285 

88.11-5749.69 

5494 

SMC X-1 

2843.66-5008.01 

704 

88.35-5747.51 

5405 

SAXJ1802.7-2017 

2698.16-5256.18 

4037 

94.24-5744.55 

2707 

XTEJ1855-026 

2704.14-5136.03 

2482 

88.37-5749.75 

5409 

Vela X-1 

2644.45-5150.95 

1665 

87.29-5747.53 

6317 

EXO 1722-363 

2668.25-5256.35 

4147 

88.11-5749.81 

4657 

OAO 1657-415 

2668.25-5256.37 

2353 

91.13-5747.53 

5065 


3. Best fit ephemerides 

3. 1. Orbital period and orbital period decay 

For each of the selected sources we first determined the mid¬ 
eclipse (superior conjunction) times, Ted, of all the eclipses 
found in the ISGRI and ASM data by using the e-fold method. 
We then improved these values by fitting together the newly de¬ 
termined epochs together with those derived from earlier obser¬ 
vations (Tn; see Appendix [A| and references therein) using the 
quadratic orbital change function; 

Tn = To+nP oi-b -I- 2 orb^orb- (1) 


Here Porb is the orbital period in days, Peak is the period deriva¬ 
tive at the epoch Tq and n is the integer number of elapsed bi¬ 
nary orbits. All the new measured orbital ephemeris are reported 
in Table We show in Fig. [T] the best fit models to the data 
obtained by assuming either a simple linear orbital evolution 
or including a quadratic orbital decay term (the residuals from 
these fits are also shown). For all those systems characterized 
by a non-negligible eccentricity (> 0.1), we fitted separately the 
mid-eclipse epochs, Ted, and the mean longitude, T„/ 2 , epochs. 
For OAO 1657-415 we fit together the epochs of Ted and 7^/2 as 
the periastron angle of the source is known to be around 90° and 
therefore Ted ~ T^r/z (see Sec. 3.2 1 . All the least square fits in this 


^ ftp://legacy.gsfc.nasa.gov/xte/data/archive/ASMP 
http://heasarc.nasa.gOv/docs/xte/recipes/asm_recipe.html#barycenter 


work were performed by using the IDF tool MPFIT ( |Markwardt| 
20091. We obtained in all cases ~ 1.0 - 1.3. 


Th the case of the five sources LMC X-4, Cen X-3, 4U 1700- 
377, SMC X-1, and 0A01657-415 we measured a significant 
orbital period derivative (see Table [^. For 4U 1538-522, the 
Porb value has been obtained excluding the mid-eclipse times 
derived from the Uhuru and Ariel lightcurves (|Cominsky & 
|M oraes 1 199 1[ [Davison et al.||1977l ); in the case of OAO 1657 


415 we excluded the orbital epoch time reported by [Barnstedt 


et al. (2008 1 . All these values were affected by relatively large 


uncertainties compared to the others and did not provide signifi¬ 
cant improvements to the fits. For the five sources characterized 
by a value of Porb consistent with zero (see Table]^, the reported 
orbital epochs and periods were determined by using the linear 
orbital change function (see Eq.[T]). 


3.2. Apsidal advance 

In eccentric orbits the time of mid-eclipse, PecU determined from 
the eclipses in the X-ray lightcurves, and the time of mean longi¬ 
tude, Tni 2 , determined through the pulse arrival times technique, 
do not coincide. In this case we distinguished the orbital period, 
Poib, 7 r/ 2 . defined as the time elapsed between two successive pas¬ 
sages at the same mean-longitude / = f, and the eclipse period, 
Rorb.eci, defined as the difference between two successive mid¬ 
eclipse epochs. The mean orbital longitude I can be expressed as 
I - M + u), where, M is the mean anomaly and to is the argument 
of periapsis. With the orbital relations, expressed in the first or¬ 
der of the eccentricity, the time delay between Ped and T „/2 can 
be written as (see e.g., ^ ~ ^ ^ ' 

references therein): 


I der Klis & Bonnet-Bidaud|1984[ and 


T, T, sPorb,M-/2 

Tjij'2 Tqq\ — cos Cl). 


( 2 ) 


For an eccentricity of e ~ 0.1 and a periastron angle of w ~ 
150° as in Vela X-1, this can result in a lag of 0.3 day. This 
periodic lag is zero for eccentric orbits if w = +90° and maximal 
when w = 0 or 180°. If the periastron angle to is constant, then 
Porb,eel = Porb,;r/ 2 - However, if the pefiastTon is moving at a rate 
io, the difference between Porb.eci and Porb, 7 r /2 at the first order of 
its eccentricity, e, can be obtained by differentiating Eq. [^ (see 
[Deeter et al.|1987] ): 


ep2 

nn orb,7rl2 . . 

-^orb,;r/2 “ -^orb,ecl 

7T 


(3) 


Equations l^andj^allow estimating to and dt if sin w 4 0. For cir¬ 
cular orbits we have Ped = T„i 2 and thus Porb.eci = Poib,n/ 2 - From 
Eq. (1^ and the value of APobs = Porb.eci “ Porb,;r /2 determined 
from me observations, we calculated also the apsidal motion to. 
The same epoch time is used to calculate both (Ped, Porb eci) and 
(P./2, PorbTi/i), as this time corresponds to the epoch of to (see 
Table[T}. ’ 

Among the four sources with e > 0.04 considered in the 
present work (see Table [^1 a detailed comparison of the mea¬ 
sured epoch time lag APobs = Pjr/2 - Ped, obtained from our 
best fit ephemerides with the predicted lag, with APcaic = 
ePorb.K/i cos to/n could not be carried out for OAO 1657-415 
(its periastron angle is w 90° and thus Ped - Tj^/i), 
XTEJ1855-026 (only one T^7ij2 ^^orb,jr/2 ^1*6 SV3.ll3.blc tO (IstC, 
see Appendix 0 and EXO 1722-363 (no reliable eccentric¬ 
ity measurement is available in the literature; Thompson et al. 


2007| ). For 4U 1538-522 we used the eccentricity, e, and the 
periastron angle, to, from |Clark| (|2000|l and [Mukheijee et al.| 
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Table 3. Updated epochs, orbital periods, and period decay of the ten sources obtained by using all the available mid-eclipse times. 
We indicated in brackets the uncertainties at Icr c.l. on the last digits of each reported value. 


Source 

To,eel (MJD) 

Foi-b.eci (days) 

Porb/Porb (10-'’yr-‘) 

oj (deg/yr) 

EMC X-4 

53013.5910(8) 

1.4083790(7) 

-1.00(5) 

- 

Cen X-3 

50506.788423(7) 

2.08704106(3) 

-1.800(1) 

- 

4U 1700-377 

53785.850(7) 

3.411581(7) 

-1.6(1) 

- 

4U 1538-522“* 

52855.061(13) 

3.7284140(76) 

-0.7(8) 

1.3(6) 

SMC X-1 

52846.688810(24) 

3.891923160(66) 

-3.541(2) 

- 

SAX J1802.7-2017 

* 52168.245(34) 

4.5697(1) 

17(29) 

- 

XTE J1855-026* 

52704.009(17) 

6.07415(8) 

-12(13) 

- 

Vela X-1* 

42611.349(13) 

8.964427(12) 

-0.1(3) 

0.41(27) 

EXO 1722-363* 

53761.695(19) 

9.74079(8) 

-21(14) 

- 

OAO 1657-415 

52674.1199(17) 

10.447355(92) 

-3.4(1) 

- 


Notes. “Epoch time and orbital period is derived from the linear fit; the orbital decay has been evaluated excluding the Uhuru and Ane/ data. * The 
orbital period derivative, P are consistent with zero and thus the epoch times and orbital periods are estimated using Eq.[^(linear fit). 


Table 4. In this table, the predicted ATcaic (Eq. |^, the measured epoch time lag, ATobs = T„i 2 - Ted, and the orbital period lag, 
APobs = Porb.eci - Poib, 7 r /2 are given for the two sources 4U 1538-522 and Vela X-1. The apsidal advance angle, io, is also reported. 
We indicated in brackets the uncertainties at Icr c.l. on the last digits of each reported value. 


Source 

Teel 

Tnjl 

P orb,ecl 


'^7'obs 

‘^7’calc 


<b 


(MJD) 

(MJD) 

(day) 

(day) 

(day) 

(day) 

(xlO^'* day) 

(deg/yr) 

4U 1538-522“ 

50450.234(11) 

50450.221(11) 

3.7284140(76) 

3.728337(22) 

-0.013(15) 

0.0905(69) 

0.72(23) 

2.3(7) 

4U 1538-522* 

52851.332(13) 

52852.3207(92) 

3.7284140(76) 

3.728382(11) 

0.011(16) 

0.164(67) 

0.32(13) 

1.3(6) 

Vela X-1 

42611.349(13) 

42611.1693(43) 

8.964427(12) 

8.9644061(64) 

-0.180(14) 

-0.2275(55) 

0.21(14) 

0.41(27) 


Notes. “Th e best fit val ues T „/2 and Povb.ii /2 were obtained including all the epochs available in t he literature until the wo rk of (the last one is 
reported by Clark 2000 1 . *Same as before but using all epochs available until the work published by Mukherjee et al. (2006 '. 


( 2006| l; for Vela X-1 we used the ephemeris and orbital param¬ 
eters from [Rappaport et al. (1976 see also Table [T]). For both 
sources the best fits are obtained using the linear orbital change 
function (Eq.[T]|. All results are reported in Tablej^ In the case of 
4U 1538-522 two different fits are reported. In the first, we fit all 
the mean longitude epochs iT„/ 2 ) until 1997 with the last one re 
ported by Clark] ((2000]l; in the second, we included all the longi 


tude epochs until 2003, with the last one reported by Mukherjee 
|et al.| ( |2006| l. The two cases give compatible results (to within 
the uncertainties) only if the epochs reported from |ClaA ( 2000| l 
are excluded from the second fit. Note that the same mid-eclipse 
time, Ted, has been used for both fits (see Table [^. 


4. Folded lightcurves 

The updated ephemerides we obtained for the ten HMXBs al¬ 
lowed folding their lightcurves with an unprecedented accu¬ 
racy. We folded for each source the dwells RXTE/ASM light 
curves in the 1.3-3 keV, 3-5 keV, and 5-12 keV bands, and 
the INTEGRAL/lSGRl light curves in the 17^0 keV and 40- 
150 keV energy bands by using 128 phase bins (see Fig. 
and[^. For EMC X-4, Cen X-3, 4U 1700-377, SMC X-1, and 
OAO 1657-415 the orbital period derivative was also taken into 
account during the folding. Given the values of PorhlPoib in 
Table we obtained for these sources a maximum variation of 
the orbital period (between the first and last mid-eclipse time 
of the RXTE/ASM dataset) of ~ 0.088, 0.158, 0.139, 0.310, and 
0.299 days, respectively. These delays have been calculated by 
using the quadratic term in Eq. [prewritten as (Af)^T’orb/T’orb (we 
replaced nPoib hy the observational elapsed time, Af). The orbital 
derivatives for these sources are not negligible, as the derived 
correction factors in time are larger than a phase bin (1 /128) and 


the inclusion of these corrections significantly improve the shape 
of the folded lightcurves. 

The fluxes (cts/s) of the ten sources measured during the oc- 
cultation of the compact object by the companion star are re¬ 
ported in Table |p In all cases, the fluxes in the lower energy 
bands (1.3-3 keV, 3-5 keV, and 5-12 keV energy bands) are con¬ 
sistent with zero. This suggests that in the soft X-ray domain the 
source emission is strongly absorbed by the extended corona of 
the companion star, as well as from the companion star itself (in 
the literature, residual fluxes have been reported at a level that is 
too low to be detectable by the ASM on-board RXTE; see, e.g., 
[Ebisawa et al.|199^|Lutovinov et al.|2000[|Vrtilek et al.|2001| l. 

At higher energies (17^0 keV and 40-150 keV), the observed 
residual X-ray fluxes are most likely due to the presence of an 
extended X-ray scattering region (e.g., the accretion disc around 
the compact object, the X-ray irradiated surface of the compan¬ 
ion star, or a cloud of material diffused around the system and 
produced by the intense wind of the supergiant star). These find¬ 
ings, together with the different shapes of the eclipse ingress and 
egress in Fig.|pand|p are discussed in Sect.|5.4| 


4.1. Semi-eclipse angles 


We used the folded lightcurves in the energy band 17^0 keV 
to estimate the semi-eclipse angle, 6^, of each eclipse observed 
from the ten selected sources in Table [p In this energy band, all 
eclipses look sharp and symmetric, thus permitting to achieve an 
unprecedented accuracy in the determination of the occultation 
time (see Fig.[^and[^. We followed the fitting method described 


Rubin et al. (|1996|l. The duration of each eclipse was calcu¬ 


lated from the measured phase of the eclipse ingress and egress. 
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Fig. 1. All measured mid-eclipses, Ted, mean longitudes, 7'^/2, and orbital epochs are shown with the best fit models as a function 
of the orbit numbers. In all cases, the best fit model is a quadratic or linear fit to the epochs. The lower panel in each figure shows 
the residual from the best fit. 


The average phase e) (defined as the phase at which 99% of 
the source flux is occulted at eclipse ingress or egress) and the 
transition phase width were estimated by fitting indepen¬ 
dently the ISGRI light curve in the 17^0 keV and (whenever 
possible) in the 40-150 keV band with the function: 


Fj = Fcons exp{ln(0.01)exp [-{(pj - 0(i,e))/T(i,e)]). (4) 


Here j is the phase bin number and Fcons is the averaged count- 
rate of the source outside the eclipse determined prior to the 
fit. The 0.5-1.0 orbital phase was considered to fit the eclipse 
ingress and the 1.0-1.5 orbital phase was used for the egress. 
The duration of each X-ray eclipse for the ten sources was cal¬ 
culated as Aifi - (1.0 - 0i) H- 0e in the (17—40 keV) band, thus 
providing the semi-eclipse angles, 0e, reported in Table (all 
uncertainties are given at Icr c.l.). We verified that the eclipse 
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Fig. 2. Folded lightcurves of the first 6 eclipsing HMXBs in TableData in the energy range 1.3-12 keV were obtained from the 
RXTE/ASM, while harder X-ray data are from INTEGRAL/ISGRl (17-150 keV). In all cases, the lightcurves have been folded by 
using the updated ephemerides obtained in the present work. 


Table 5. Measured cts/s during the occultation phases of the ten 
selected sources in this work. Note that 0.2 cts/s (0.1 cts/s) in the 
17^0 keV (40-150 keV) energy band in ISGRl correspond to 
roughly 1 mCrab. We indicated in brackets the uncertainties at 
90% c.l. on the last digits of each reported value. 


Source 

17-40keV 

cts/s 

40- 150keV 
cts/s 

LMC X-4 

0.9(1) 

0.7(1) 

Cen X-3 

0.27(8) 

- 

4U 1700-377 

1.31(8) 

0.50(6) 

4U 1538-522 

- 

- 

SMC X-1 

0.48(16) 

- 

SAX J1802.7-2017 

0.35(7) 

- 

XTEJ1855-026 

- 

- 

Vela X-1 

1.02(7) 

0.15(6) 

EXO 1722-363 

- 

- 

OAO 1657-415 

0.60(7) 

0.23(7) 


All the semi-eclipse angles we measured and reported in 
Table are a few degrees shorter than the values reported in 
the literature, the only two exceptions being 4U 1700-377 and 
LMC X-4. In the case of 4U 1700-377, our measured ingress 
phase is shorter than that measured by Rubin et al. (1996 1 , thus 
translating into a slightly longer eclipsing phase. Eor LMC X-4 
our determined semi-eclipse angle of Og ^ 16° is around half 
the value ~ 27° reported in the literature ( |Li et al.|1978 1. This is 
most likely due to the very few number of eclipses available in 
the past and the different energy band used to estimate 6^ (| White 


1978t Pietsch et al. 1985|l. The usage of the long term obser¬ 


vations available at present allow reducing the temporal variable 
distortions of the eclipse profile and provide a more reliable mea¬ 
surement of the semi-eclipse angle and eclipse duration. 


5. Discussion 

5.1. Masses of the 10 Neutron Stars 


Determining the equation of state (EoS) of matter at densities 
comparable to those inside NSs is one of the most challenging 
semi-angles measured in the 17^0 keV and 40-150 keV energy problems of the modern physics and can only be addressed based 
band are compatible to within the uncertainties. on observations of astrophysical sources. Models proposed in the 
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Fig. 3. Same as Fig. I^but for the remaining 4 sources in 



past years can be tested against observational results especially 
evaluating the maximum NS mass that each EoS model is able 
to sustain (see e.g., Lattimer & Praka^|2001 1 . Very soft EoSs 
predict maximum NS masses in the 1.4—1.5 Mq range (this oc¬ 
curs when the NS core is made by exotic matter as kaons, hyper- 
ons, and pions), whereas stiff EoSs can reach up to 2.4—2.5 Mq. 
More massive NSs can thus provide stronger constraints on the 
EoS models. As discussed by Rappaport & Joss (19831, eclips¬ 
ing HMXBs hosting X-ray pulsars provide a means to measure 
the NS mass and thus put constraints on their EoS. 


In eclipsing HMXB the parameters needed to solve the 
equations that lead to the determination of the NS mass are: 
the system orbital period Forb. the projected semi-major axis 
a^sini, the eccentricity e, the periastron angle a>, and the du¬ 
ration of the eclipse expressed as the semi-eclipse angle 6^ 
(see Sec. |4. l| i. The semi-amplitude of the NS radial velocity 
can be inferred from the results of timing analysis providing 
Xx = Ina^sini /(Porhi^ - The semi-amplitude of the 

radial velocity of the optical component, Xopt, can be deter¬ 
mined by optical and/or UV spectra of this star (see Table 
we note that the values of Xopt determined through these tech¬ 
niques might be affected by uncertainties related to the compan¬ 


ion st ar modelling see, e.g., Hammerschlag-Hensberge et al. 
2003 Abubekerov| 2004[ Koenigsberger et al.||2012| l. The NS 


radial velocities, Xx, reported in Table |7| are inferred using the 
values of a^sini, e, and Port reported in Table and Two ad¬ 
ditional key parameters are required and can be estimated from 
theoretical arguments: (/) the Roche lobe filling factor p (i.e., 
the ratio of the supergiant’s radius to that of its Roche lobe 
P - Ro-ptIRi) and (//) the ratio of the spin period of the super¬ 
giant to its orbital period Q. 

Eor sake of completeness, we review below the main equa¬ 
tions to determine the NS masses (see also [Rappaport & Joss 


19831. The masses of the optical supergiant companion, Mopt, 
and the neutron star, M^, can be written in terms of the mass 
functions: 


34opt — 


x3p„,b(l-ey/2 

In G sin^ i 


(l+qf 


(5) 


^ The strong stellar wind, which is responsible for the mass transfer, 
makes the task of generating a radial velocity curve of the donor star 
difficult. Another issue that can result in difficulties is also when the 
compact companion heats up one side of the donor star, thereby distort¬ 
ing spectral line profiles, i.e.. X-ray heating take place. 
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and 




^opt^orbd 


2n G sin^ i 


,2)3/2 

(1 + 




(6) 


where the mass ratio q is defined as, q = M^/Mopt = KgptlK^. 
Assuming a spherical companion star, the inclination angle of 
the system, i, is found to be related to the semi-eclipse angle 
through the geometric relation: 


Sinz - 

cos Oq 


(7) 


We also made use here of the equation p - /?opt/^L- The ratio 
between the Roche lobe radius and the separation of the two 
stars, can be approximated as Poss & Rappaport|1984| l: 

— ^ A + B log q + C log^ q, (8) 

a 

where: 

A = 0.398 - 0.026Q^ + 0.004Q^ 

B = -0.264-H 0.052^2 - O.OISQ^ 

C = -0.023 - 0.005GI (9) 


For HMXBs with circular orbits the ratio Ri^ja is constant. 
When the orbit is eccentric the Roche Lobe filling factor, j0, is de¬ 
fined at the periastron, R]^, and the separation between the cen¬ 
ters of masses of the two stars varies with the orbital phase, n 
these cases, the separation between the centers of the two stars 
at mid-eclipse time is given by a' - a{l - e^)/(l -H ecosu), 
where oj is the argument of periastron reported in Table [T] In 
our calculation the Roche lobe radius is estimated at the eclipse 
phase. The approximated Roche lobe radius, Ry^, is determined 
with an accuracy of about 2% over the ranges of 0 ^ G ^ 2 
and 0.02 ^ ^ ^ 1 (Joss & Rappaport 1984 1 . Once the full 
set of input parameters above are given, the values of i, M^, 
Mopt, a, R]^, and Ropt can be determined at Icr c.l. by means of 
Monte-Carlo (MC) simulations. Uncertainties on all parameters 
are evaluated in the simulations by assuming that their values 
follow a Gaussian distribution. We used a linear distribution for 
the Roche lobe filling factor,/!, spanning the range 0.9-1.0. Note 
that for all the HMXBs considered here it is known from their 
optical lightcurves that the supergiant stars hosted in these sys¬ 
tems are nearly filling their Roche lobes (see e.g., |Tjemkes et al.| 
|1986| and references therein). 

Following the technique described above, we first computed 
for the four sources with known Q (LMC X-4, Cen X-3, SMC 
X-1, and Vela X-1) the remaining output parameters reported 
in Table In order to verify the reliability of this method, we 
compared the observational reported projected stellar rotation 
velocity, Vrotsin i, measurements with our calculated value using 
results in Table and the following equation: 


Ikms-WldaysA/^o/ (l+e)i/2- 


( 10 ) 


In Table 1^ we compare the calculated and measured Vrotsini val¬ 
ues. These values agree well within the uncertainties. Based on 
these results, we assumed for the sources 4U 1538-522, SAX 
J1802.7-2017,XTEJ1855-026,EXO 1722-363, andOAO 1657- 
415 a mean value of G = 0.91(20). As these systems are not 
circularized, Q. is unlikely to be too close to unity. If tides are 
efficient, then we might expect that the rotation of the donor is 


Table 8. Comparison between the calculated and measured val¬ 
ues of Vrotsin i. 


Source 

ViofSin / 

km 

Vrotsin r 

km 

(JS = 0.9- 1.0) 

LMC X-4 

240(25)'' 

257.9(7) 

Cen X-3 

200(40)“ 

204.8(8) 

4U 1538-522 

180(30)* 

224.9(7) 

SMC X-1 

170(30)“ 

178.0(8) 

SMCX-1 

172(1.5)“ 

178.0(8) 

Vela X-1 

116(6)'* 

130.3(2) 


Notes. 


der Meer et al.|l^07t, ^va n Kerkwijk ^tj^| l995a| see 


reference therein), jReynolds et al. ( 1993j , ‘ Zuiderwij^ ( |l995| . *This 
work, using Eq.|10|and our results reported in Table |7| and Table|^ 


synchronized at periastron, where the tides are most effective 
( Hut|198r| ). In this case, the star’s angular velocity is larger than 
the orbital velocity at mid-eclipse, thus leading to Q > 1. Eor 
this reason we increased the uncertainty on Q to include values 
slightly larger than unity. 

4U 1700-377 was considered separately, as its projected 
semi-major axis OxSin i is unknown and thus the corresponding 
Xx cannot be determined as we have done for all other sources. 
Eor this source, we first estimated Q from Eq. [T^by using the 
observationally measured value of Vitsin/, Rorb. and the com¬ 
panion star radius Ropt = 21.9(1.9)Ro (see Table |7| and Clark 
et al.|200^ . We thus performed several MC simulations with a 
variable Xx, until the Ropt obtained from one of the simulations 
was compatible (to within the uncertainties) with the observa¬ 
tional value. The outputs of this simulation were then used to 
fill all other relevant parameters for the source 4U 1700-377 in 
Table |2] 

Our findings on all neutron star masses for the considered 
HMXBs are summarized in Fig.|^ together with previously pub¬ 
lished values. The neutron star masses for LMC X-4, Cen X- 
3, 4U 1538-52, SMCX-1, and Vela X-1 are from [Rawls elaL] 
( |201 l| l, while the mass of the neutron star hosted in 4U 1700-377 
is derived from|Clark et al.|(|20021l. The corresponding values for 
SAXJ1802. 7-2017,EXO 1722-363, and OAO 1657-415 are ob¬ 
tained from [Mason et al.| ( |2011[|2010| |2012| l. We note that our 
estimated uncertainties on the neutron star masses are somewhat 
more conservative than those reported b y[Rawls et al.[(|201 l| l and 
slightly smaller than those obtained by Clark et al. ( 2002| i and 
[Mason et al.[ ( [201 l1[2010|[2012| l. The reason behind these differ¬ 
ences is that these authors derived their uncertainties either from 
analytical calculations or through numerical simulations includ¬ 
ing a particularly refined treatment for the Roche-Lobe size of 
the supergiant star. In the present work, all uncertainties are, 
instead, derived directly from the MC simulations and account 
for the most recently updated system parameters obtained from 
multi-wavelength observations. In particular, we reported for all 
sources in Table the most accurately measured values of the 
semi-eclipse angles, which play a crucial role in the MC simula¬ 
tions. We note that in all cases our measured semi-eclipse angle 
is smaller than reported values in the literature (see Table|^, and 
thus the NS masses we estimated are generally larger. This seems 
reasonable especially for4U 1538-522 and SMCX-1, as the NS 
masses estimated before where <1 Mq. The only exception to 
this trend is 4U 1700-377, for which our estimated semi-eclipse 
angle is comparable to that obtained by [Rubin et al.[(fT996ll. 


8 

































































Falanga et al.: The ephemeris, orbital decay, and masses of 10 HMXBs 


Table 6. Best fit orbital profile parameters in units of the orbital phase. For each source we indicated with 6^ the semi-eclipse angle 
measured from our ISGRI lightcurves (17^0 keV energy band), e is the dehned as the phase at which 99% of the source flux 
is occulted at eclipse ingress or egress, and Ti e the corresponding transition phase width. The values of the semi-eclipse angles 
reported in the last column, 0e,oid, are from the literature. We indicated in brackets the uncertainties at Icr c.l. on the last digits of 
each reported value. 


Source 

0. 

Ti 

<Pe 

Te 

(deg) 

^e,old 

(deg) 

LMC X-4 

0.951(4) 

0.020(1) 

0.039(3) 

0.012(2) 

15.8(8) 

27(2)" 

Cen X-3 

0.922(1) 

0.019(1) 

0.077(1) 

0.0187(6) 

27.9(3) 

33(1)* 

4U 1700-377 

0.912(4) 

0.0185(5) 

0.0894(4) 

0.0126(2) 

32(1) 

29(2)" 

4U 1538-522 

0.936(4) 

0.018(2) 

0.052(4) 

0.0147(5) 

21(1) 

29(2)'' 

SMC X-1 

0.939(9) 

0.0139(4) 

0.066(5) 

0.007(2) 

23(2) 

28(2)" 

SAX J 1802.7-2017 

0.94(1) 

0.038(9) 

0.11(1) 

0.012(7) 

31(2) 

35(5)" 

XTE J1855-026 

0.920(5) 

0.0129(3) 

0.097(6) 

0.016(4) 

32(1) 

42(6)' 

Vela X-1 

0.9202(4) 

0.0209(1) 

0.0899(2) 

0.0084(2) 

30.5(1) 

33(3)s 

EXO 1722-363 

0.94(2) 

0.013(8) 

0.088(4) 

0.004(2) 

26(4) 

32(2)* 

OAO 1657-415 

0.933(1) 

0.0184(6) 

0.049(2) 

0.0352(1) 

20.9(4) 

30(1)' 


Notes, ‘van der Meer et al.l(l2007l and references therein), IClark et al.l(ll988i. iRubin et al. 1 1996b, 

van Kerkwiik et al. 

1995a 1 , 

'Hill et al. 

2005 1 , ^Corbet & Mukai ( 2002^, jvan Kerkwijk et al. 

1995bb, 'ICorbet et al. ( 2005b, Chakrabarty et al. 

119931. 


Table 7. Input and output parameters for the MC simulations used to finally estimate the mass of the NSs hosted in the 10 HMXBs 
considered in this work. We indicated in brackets the uncertainties at Icr c.l. on the last digits of each reported value. 


Source 


Input 





Output 




K, 

Kopt 

9 

n 

i 

Mn, 

^opt 

^opt 

Rh/a 

a 

km 

km s ^ 



deg 

Mo 

Mo 

Re 


Re 

LMC X-4 

407.8(3) 

35(2)" 

0.086(4) 

0.97(13)" 

59.3(9) 

1.57(11) 

18(1) 

7.4(4) 

0.59(1) 

14.2(2)' 

Cen X-3 

414.317(9) 

28(2)" 

0.066(5) 

0.75(13)" 

65(1) 

1.57(16) 

24(1) 

11.4(7) 

0.63(1) 

20.2(4) 

4U 1700-377" 

435(10) 

19(1)* 

0.043(2) 

0.47(4)"" 

62(1) 

1.96(19) 

46(5) 

22(2) 

0.694(6) 

35(1) 

4U 1538-522 

316(10) 

20(3)" 

0.06(1) 

0.91(20)'" 

67(1) 

1.02(17) 

16(2) 

13(1) 

0.53(3) 

22(1)' 

SMC X-1 

299.631(5) 

20(1)" 

0.067(4) 

0.91(20)" 

62(2) 

1.21(12) 

18(2) 

15(1) 

0.61(2) 

21.9(1)‘ 

SAXJ1802.7-2017 

324(5) 

24(3)'' 

0.07(1) 

0.91(20)" 

72(2) 

1.57(25) 

22(2) 

18(1) 

0.61(2) 

33(1) 

XTE J18f55-026 

289(5) 

20(3)" 

0.07(1) 

0.91(20)" 

71(2) 

1.41(24) 

21(2) 

22(2) 

0.63(3) 

40(1)' 

Vela X-1 

278.1(3) 

23(2)" 

0.081(5) 

0.67(4)*-' 

72.8(4) 

2.12(16) 

26(1) 

29(1) 

0.595(6) 

59.6(7)' 

EXO 1722-363 

226(2) 

25(5)^ 

0.11(2) 

0.91(20)" 

68(2) 

1.91(45) 

18(2) 

26(2) 

0.58(3) 

52(2) 

OAO 1657-415 

222.60(6) 

21(4)« 

0.09(1) 

0.91(20)" 

67.9(9) 

1.74(30) 

17.5(8) 

25(2) 

0.52(2) 

53.1(8)' 


Notes. 


van der Meer et al. 


,lHammerschlag-Hensberge et al. i 2003 1 ;' van Kerkwijk et al. jl995a see reference therein). ‘^'^Vmtsin i = 150( 10) 


km/s (Clark et al, 2002); ^Mason et al.| ( |2011[ >; ' iQuaintreil et al.| i 2003i; ^Mason et al.| l 20101, 'Mason et al. j2012^ ; ' Howarth et al.| 1 1997 1 ; 


Zuiderwijk| i 1995 i ;' evaluated at periastron; '"we used the same fl as for SMC X-1 since the values of Vrotsin i for these systems are compatible 
to within the uncertainties (see also Rawls et al.|201 l|l; "we assumed a value of derived from the average of the other systems in this table 
with properties close to XTE J1855-026; “we assume ai conservative fl value in order to be close to the Roche lobe radius; * we treat 4U 1700-377 
differently, since the project semi-major axis, a^sin i, is unknown; see Sec. mi 


5.2. Orbital period change 

Five out of the 10 HMXBs analyzed in the present work showed 
evidence for a signihcant orbital period decay, i.e., LMCX-4, 
CenX-3, 4U 1700-377, SMCX-1, and OAO 1657-415 (our val¬ 
ues can be compared with previous measurements published in 
the literature by using the references reported in the Tables in 
Appendix [All. For all these sources, a number of different mech¬ 
anisms to explain the orbital period decay has been invoked 


(see e.g., |Kelley et al.||1983| van der Klis 

983 van der Klis 

& Bonnet-Bidaud||1984 

ILevine et al.||1993 

Rubin et al.||19961 

Safi Harb et al. 19961 1] 

^evine et al. 2000 

enke et al.| 2011b. 


The hypothesis that the period change is due entirely to angu¬ 
lar momentum loss in the stellar wind, i.e. mass transfer with 
no interaction, can be ruled out. More realistic models proposed 
that the orbital decay is mostly driven by the tidal interaction 
and the rapid mass transfer between the two objects. The latter 
takes place due to the fast winds that characterize most of the 


donor stars in HMXBs. In these models, the asynchronism be¬ 
tween the orbit and the rotation of the donor star is maintained 


by the evolutionary expansion of the donor star (see e.g., Lecar 
al.|1976HSafi Harb etal.| 199^ [Levine et al.|2000[|van der Klis 


& Bonnet-Bidaud 1984[ and references therein). If the donor star 
is rotating with an angular velocity close to synchronization, the 
matter ejected leaves the system and carries an extra momentum 
due to this rotation. This loss of rotational momentum affects 
the orbital period only if there is a coupling between the rotation 
of the donor star and the orbit of the neutron star. Tidal interac¬ 
tions could synchronize the orbit with the rotation of the donor 
star and allow the transfer of rotational momentum of the pri¬ 
mary star to the orbital momentum. Tidal effects on the radial 
velocity curve of the donor star in Vela X-1 have been investi¬ 
gated deeply, but not final conclusions could be drown due to 
the limited quality of the data ( Koenigsberger et al.|2012 1 . Tidal 
interactions in close binary systems has heen investigated in de¬ 
tail also theoretically by using different approaches, e.g., assum- 
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LMC X-4 
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I 

I 
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Vela X-1 
EXO 1722-363 
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Fig. 4. The masses of the ten eclipsing HMXBs. The neutron 
star masses determined in this work are shown with solid lines. 
Values from the literature are represented with dashed lines. The 
error bars correspond to uncertainties at Icr c.l. The dashed ver¬ 
tical line indicate the canonical neutron star mass of 1.4 Mq. 


ing the rotation of the donor star is synchronized or pseudosyn- 
chronizecl^with the orbital motion of the compact star (see e.g., 
[Moreno et al.T20TTl[HutnMTl|Zahn|1977||Lecar et al.|1976| l. A 

number of other works argued that the combined effects of the 
evolutionary growth of the moment of inertia (expansion of the 
donor star) and mass loss by stellar wind from a donor star can 
explain the orbital decay in HMXBs ( Bagot|1996[ Levine et al. 
1993^. In these models, it is considered that, as the companion 


star evolves, it expands and increases its moment of inertia, lead¬ 
ing to a slow down of the star rotation. Tidal torques would then 
transfer orbital momentum to synchronize the binary system, fi¬ 
nally causing an orbital decay (see e.g., Levine et al.|1993 and 
references therein). As the lifetime of a supergiant star does not 
exceed a few 10^ yr, we expect that most of the systems consid¬ 
ered in this work are nearly synchronized but not fully circular¬ 
ized yet. This is in agreement with the non-negligible values of 
the eccentricity measured in some of the eclipsing HMXBs (see 
Tableg. 

The updated and newly reported estimates of the orbital pe¬ 
riods decay of the ten HMXBs considered in the present work 
can be used as inputs in future works on these topics. We expect 
that this might help discriminating among the different mecha¬ 
nisms proposed to explain the orbital period decay in HMXBs 
and to estimate their orbital circularization and synchronization 
timescales, as well as to discuss the Darwin instability (see e.g., 
Lai et al.|1994|l for these systems. Additional calculations on all 


these effects are beyond the scope of the present paper. 

In order to help future works to evaluate numerically differ¬ 
ent theoretical models, we finally report for reference also the 
stellar wind velocities and mass loss rates of all sources con- 


® Pseudosynchronized systems are binaries in which the donor star is 
synchronized only at the periastron, where the tides are more effective 


(Hut 19811. 


Table 9. Properties of the stellar winds in the 10 HMXBs consid¬ 
ered in this paper as inferred from observations and theoretical 
calculations. 


Source 

Vco 

km 

Observed 

Vco 

km s~^ 
Theoretical* 

m 

Mo yr-‘ 
Theoretical* 

EMC X-4 

1350(35)' 

1950(600) 

2.4 X 10"' 

Cen X-3 

-500^ 

2050(600) 

5.3 X 10-3 

4U 1700-377 

1700(100)3 

1850(550) 

> 2.1 X 10-® 

4U 1538-522 

_4 

1000(300) 

8.3 X 10-3 

SMC X-1 

_4 

870(260) 

1.5 X 10-'’ 

SAX J1802.7-2017 

4 

680(200) 

6.3 X 10-3 

XTE J1855-026 

_4 

620(190) 

<(0.2- 1.1) X 10-3 

Vela X-1 

600(100)3 

640(190) 

< (1.0-5.3) X 10-3 

EXO J1722-363 

_4 

650(200) 

9.0 X 10-3 

OAO 1657-415 

~2503 

200(60) 

< (1.1 -5.6) X 10-3 


Notes. References: * [Boroson et al. 1 19^; ^ Wojdowski et al. 1 2003 i; 


van Loon et al.j l 2001 1 ; No consolidated measurement availa ble m the 
literature - at the best of our knowledge; ^ 


Mason et al. 


20121 . 


sidered in this paper in Table These have been calculated by 
assuming the spectral types reported in Table [T] In Table we 
indicate both the terminal wind velocities inferred from observa¬ 
tions in the literature and those calculated theoretically. For the 
latter we followed [Martins et al. (2005 and [Kudritzki & Puls 


( 2000 ( for O-stars, an d [Markova & Puls ( [2008| l, Crowther et al. 


( 2006 (, [Lefever et al.[ ( [2007| l for B-stars. In all cases, the nom 


inal mass loss rates have been reduced by a factor of ~ 3 in 
order to account for clumping and uncertainties on the terminal 
wind velocities are assumed to be at the 30% ( [Markova & Puls[ 
2008[ Repolust et al. 2004| l. As all data in Table have been 


calculated assuming the wind is emitted by an isolated massive 
star, it is possible that in a number of systems the disturbance of 
the neutron star lead to substantially different wind properties. 
Detailed numerical simulations have showed that this could be 
the case especially for short orbital period systems with high X- 
ray luminosities ( [Watanabe et al.|2006[ Manousakis et al.|2014| l. 
Unfortunately, direct measurements of the wind properties for 
the systems considered here are challenged by the relatively high 
absorption local to the sources and their large distance ([Chaty[ 
[20131 ). 


5.3. Apsidal motion 


Beside asteroseismology (see, e.g., Dupret et al.|2011] and ref¬ 
erences therein), the measurement of the apsidal advance in ec¬ 
centric binary systems offer an alternative possibility to investi¬ 
gate the internal structure of stars compared to theoretical cal¬ 
culations and numerical modelling. The rate at which the longi¬ 
tude of periastron of an eccentric orbit advances can be related 
to the absidal motion constant, which depends on the model of 
the mass distribution inside the star ( [Kopal|1978)|Batten|1973| l. 
X-ray pulsars in HMXBs have been far considered the most 
promising laboratories to measure accurately the apsidal ad¬ 
vance. The neutron star in these systems can be well approxi¬ 
mated as a point source with a negligible absidal motion con¬ 
stant compared to that of the massive companion. The latter is 
thus dominating the apsidal advance ( [Rappaport et al.|1980| l. 

The long term analysis carried out in this paper, allowed us 
to provide the most accurate apsidal advance measurement for 
VelaX-1 and 4U 1538-522. (note, however, that for the former 
source the measurement is only marginally significant at 1.5cr). 
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In the case of VelaX-1, the theoretical estimated apsidal mo¬ 
tion is u)~ 0.4° yr“', assuming a companion star surface temper¬ 
ature of 25000 K (as predicted accordingly to its spectral type; 
see, e.g., |Avni & Bahcall|1975[|Conti|1978l l. The measurement 
reported in Table 14 is thus fully consistent with the theoretical 
expected value (see Deeter et al.|1987 and references therein). 
To the best of our knowledge, this is the hrst time that such match 
is reported, given the reduced uncertainty on the measured value 
of (b (even though this is still only marginally signihcant). 

In the case of 4U 1538-522, no hrm measurement of the apsi¬ 
dal motion was reported in the literature. We notice that the ten¬ 
tative value provided by Mukherjee et al. (2006i was estimated 


by comparing only two measurements of cu collected in 1997 
and 2003 and thus is not comparable to the rehned estimate re¬ 
ported in the present work. Our measurement thus provides the 
hrst consolidated estimate of u for 4U 1538-522. 


5.4. Eclipse asymmetry 

The updated ephemeris for all the 10 HMXBs in Tablej^that we 
obtained in this work allowed us to fold with an unprecedented 
accuracy the long term ASM and ISGRI lightcurves of these ob¬ 
jects in different energy bands. The results are shown in Fig. 
and The available broad-band energy coverage permitted to 
study in detail the averaged prohle of the eclipse ingress at egress 
for each object. We note that in the case of SAX J1802.7-2017, 
XTE J1855-026, andEXO 1722-363, this study could not be eas¬ 
ily carried out at low energies due to the strong absorption affect¬ 
ing the X-ray emission from these sources (this is expected ac¬ 
cording to their classihcation as “highly obscured HMXBs”; see, 
e.g. Chaty|20r3 for a recent review). In most of the other sources 
we found peculiar asymmetries that are more pronounced at the 
lower energies and slightly more enhanced in objects character¬ 
ized by a non-negligible eccentricity. Remarkable asymmetries 
are recorded for the eclipse ingresses and egresses of 4U 1700- 
377 and VelaX-1. Evident eclipse asymmetries are also dis¬ 
played by 4U 1538-522 and OAO 1657-415. 

In most of the sources considered in this work, single eclipse 
prohles have been studied in the past by using focusing X-ray in¬ 
struments capable of high resolution spectroscopy (e.g., XMM- 
Newton and Chandra; see [Watanabe et al.|2006i and references 
therein). These observations proved to be particularly powerful 
to analyze how the X-ray emission from HMXBs is affected 
by the presence of inhomogeneities in the stellar winds (|Sako 


et al. 


2003|l. Such inhomogeneities (“clumps”, see e.g.|Puls et al. 


2008 for a recent review) are transported away from the star 
with typical velocities of a smooth wind and thus can induce 
changes in the mass accretion rate (and X-ray luminosity) on 
timescales as short as 100-1000 s. The feedback of X-ray radia¬ 
tion is also known to affect the structure of the wind on similar 
time scales, as high energy photons ionize the metal ions in the 
wind and dramatically reduce the wind acceleration mechanism 
(in turns affecting the velocity and density profile of outflows 
from the massive stars). 

The eclipse prohles we discussed above provide, instead, im¬ 
portant information on the interaction between the compact ob¬ 
ject and the stellar wind on much longer time scales (~yrs). The 
inhuence of wind inhomogeneities is in all cases averaged away 
by the long-term integration, and only the effect of X-ray irra¬ 
diation on the wind averaged on hundreds of orbital revolutions 
can be observed from these curves. Under these assumption, we 
argue that the eclipse asymmetries are most likely due to the 
presence of “accretion wakes”. These structures are commonly 
observed in HMXBs and form as a consequence of the compact 


object intense gravity and conspicuous X-ray emission (see, e.g., 
Eig. 5 in|Blondin et al.|(|1990|l and Eig. 1-5 in|Manousakis et al. 
( I2&T2I 1). 

As the accretion wake usually trails the neutron star during 
its orbital revolution around the companion and can lead to an 
increase of the wind density around the compact object by a fac¬ 
tor of ~10-100 ( Manousakis et al.|2014 i, it is expected that ab¬ 
sorption column density in the direction of the X-ray source pro¬ 
gressively increases before the occurrence of the X-ray eclipse. 
During the egress from the eclipse, the accretion wake is located 
beyond the compact object along the line of sight to the observer 
and thus do not lead to any apparent increase of the local absorp¬ 
tion column density. 

An interesting possibility is also that this effect is enhanced 
in systems endowed with a non-negligible orbital eccentricity 
due to the different degree of ionization of the stellar wind mate¬ 
rial when the neutron star approaches or recedes from the com¬ 
panion. It is well known that in case of eccentric orbits, the ab¬ 
solute value of the vectorial sum of the neutron star orbital ve¬ 
locity and the stellar wind velocity reaches a minimum while the 
compact object approaches the companion, and then increases 
toward the upper conjunction (see, e.g., Pucci et^|2010 and 
references therein). In a wind accreting system, the amount of 
material that a neutron star can capture from the companion at 
any time is roughly comparable to where 


Ra 


2GMns 


( 11 ) 


rel 


and Vrei is the relative velocity mentioned above. It is thus ex¬ 
pected that a larger X-ray luminosity is released at lower values 
of Vrei due to the enhanced accretion rate (see, e.g., Bozzo et al. 
|2008[ and references therein). As the ionization of the wind is 
proportional to the X-ray flux and acts to substantially inhibit 
the wind acceleration mechanism, a larger density would also be 
expected around a bright X-ray emitting neutron star approach¬ 
ing the companion. Eurthermore, the X-ray flux could be evap¬ 
orating additional material from the massive star around this or¬ 
bital phase, possibly enhancing the absorption column density 
in the direction of the source during the eclipse ingress (Blondin 
& Woo|1995 1. We note that in short orbital period systems, the 


formation of (at least) temporary accretion disks when the neu¬ 
tron star is at the closest distance from the companion could also 
provide enhanced accretion rates and thus positively contribute 


to increase the mass density around the neutron star (e.g., Ducci 


et al. 20101. A quantitative estimate of all these effects would 


require more detailed calculations combining also the effect of 
the neutron star spin and magnetic held. 


6. Conclusions 

In this work, we performed orbital mid-eclipse time measure¬ 
ments on the ten known eclipsing HMXBs using publicly avail¬ 
able pre-processed RXTE/ASM and INTEGRAL/lSGRl data. Eor 
each source, we determined several new orbital mid-eclipse 
times and used them together with all previously published value 
(to the best of our knowledge) to obtain the most updated and ho¬ 
mogeneous set of ephemerides to date. The latter allowed us in 
particular to update the orbital periods of all the 10 binaries, as 
well as the decay of these periods. In hve sources (LMCX-4, 
CenX-3, 4U 1700-377, SMCX-1, and OAO 1657-415) we mea¬ 
sured a signihcant orbital period decay, with values in the range 
1.0 - 3.5 X 10“® yr“*; in all other case no signihcant decay was 
revealed even though more than 30 yrs of X-ray observations 
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were used. For the eccentric systems VelaX-1 and 4U 1538-522 
we accurately determined for the hrst time a signihcant apsidal 
advance, comparing it with previously expected values. 

Using our best estimated mid-eclipse epoch, orbital period, 
and period derivative for each source, we folded the lightcurves 
of the ten binaries obtained from the long-term monitoring car¬ 
ried out with the ASM on-board RXTE (1.3-12 keV) and ISGRI 
on-board INTEGRAL (17-150 keV). We discussed the asym¬ 
metric prohles of the eclipse ingresses and egresses, comparing 
them with theoretical expectations from simulations of accreting 
wind-fed systems. 

Einally, we used all above mentioned results within the 
MC simulations to consistently and systematically evaluate the 
masses of the NS hosted in the 10 eclipsing HMXBs.The results 
reported in this paper constitute an important database for popu¬ 
lation and evolutionary studies of HMXBs, as well as theoretical 
modelling of long-term accretion in wind-fed X-ray binaries. 
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Appendix A: Ephemeris tabies 


LMC X-4 


Epoch Time 

Satellite 

Reference 

(MJD) 




42829.494(19) 

SAS-3 

Kelley et al 

U19831 

45651.917(15) 

EXOSAT 

Uennerl 

T9 



45656.1453(8) 

EXOSAT 

Dennerl 

vm 


46447.668(11) 

EXOSAT 

Uennerl 

1991 


46481.467(3) 

EXOSAT 

Uennerl 

1991 


47229.3313(4) 

Ginga 

Woo et a 


m 


47741.9904(2) 

Ginga 

Levine et a 

Tfrai 

48558.8598(13) 

ROSAT 

Woo et al. 



51110.86579(20) 

52648.813(14) 

53013.590(14) 

53016.40(2) 

54262.825(8) 


RXTE 

INTEGRAL 

INTEGRAL 

INTEGRAL 

RXTEjhSM 


Levine et al.|j2000 
Present Work 

Present Work 
Present Work 
Present Work 


Cen X-3 


Epoch Time 


Satellite 


Reference 


40958.34643(45) 

Uhuru 

Fabbiano & Schreier 

1977 

41077.31497(15) 

Uhuru 

t^abbiano & iichreier 

wn 

41131.58181(29) 

Uhuru 

fabbiano & Sctireier 

1977 

41148.28051(16) 

Uhuru 

Fabbiano & Schreier 

1977 

41304.81533(14) 

Uhuru 

Fabbiano & Schreier 

wn 

41528.1401(3) 

Uhuru 

Fabbiano & Schreier 

wn 

41551.09798(17) 

Uhuru 

Fabbiano & Schreier 

wn 

41569.88199(11) 

Uhuru 

Fabbiano & Schreier 

wn 

41574.05610(13) 

Uhuru 

Fabbiano & Schreier 

wn 

41576.1433(1) 

Uhuru 

Fabbiano & Schreier 

wn 

41578.23037(7) 

Uhuru 

Fabbiano & Schreier 

wn 

41580.31722(9) 

Uhuru 

Fabbiano & Schreier 

1977 

41584.49193(10) 

Uhuru 

Fabbiano & Schreier 

1977 

41590.75328(15) 

Uhuru 

Fabbiano & Schreier 

1977 

41592.84025(15) 

Uhuru 

Fabbiano & Schreier 

1977 

41599.10212(15) 

Uhuru 

Fabbiano & Schreier 

1977 

41601.18930(14) 

Uhuru 

Fabbiano & Schreier 

1977 

41603.27671(21) 

Uhuru 

Fabbiano & Schreier 

1977 


42438.128(3) 

Ariel-V 

luohviiy/bi 

42786.6755(7) 

Cos-B 

van der RIis et al. 

SO 

OO 

o 

43112.26642(40) 

SAS-3 

Kellev et al. 

T9S3 


43700.83275(43) 

HEAO-1 

Flowe et al. 

1983 


43869.88910(2) 

SAS-3 

Kellev et al. 

1983 


44685.94760(5) 

Hakucho 

Murakami et al. 1 

J83 

45049.1025(1) 

Hakucho 

INagase et al. 



45428.95421(5) 

Tenma 

JNagase et al. 

1984 


47607.8688(8) 

Ginga 

JNagase et al. 

1992 



50506.788423(7) RXTE 


Kaicnur & Paul ' 2D10 


50782.279(8) 

RXTE/ASM 

Present Work 

52180.589(8) 

RXTE/ASM 

Present Work 

53136.455(14) 

INTEGRAL 

Present Work 

53574.711(8) 

RXTE/ASM 

Present Work 

54144.471(14) 

INTEGRAL 

Present Work 

54966.745(14) 

RXTE/ASM 

Present Work 

4U 1700-377 

Epoch Time 

Tnjl 

Satellite 

Reference 


41452.64(1) 

Uhuru 

Jones et al.Ul973l 


42609.25(1) 

Copernicus 

branduardi et al. 

11978 


42612.646(10) 

Copernicus 

branduardi et al. 

1978 


43001.604(10) 

Copernicus 

branduardi et al. 

1978 


43005.00(1) 

Copernicus 

branduardi et al. 

1978 


46160.840(3) 

EXOSAT 

Haberl et al. 




48722.94(31) 

Granat 

Kubm et al. 

rwB 


48651.365(31) 

BATSE 

Kubin et al. 

1996 

■ 

49149.425(27) 

BATSE 

Kubm et al. 

1996 

■ 

50783.650(15) 

RXTT’/ASM 

Present Wor 





52175.579(15) 

52861.29(2) 

53270.69(2) 

53574.341(15) 

53472.00(2) 

53785.82(3) 

54164.55(2) 

54341.97(3) 

54962.83(2) 


RXTEjASM 

INTEGRAL 

INTEGRAL 

fiXTE/ASM 

INTEGRAL 

INTEGRAL 

INTEGRAL 

INTEGRAL 

RXTEjASM 


Present Work 
Present work 
Present work 
Present Work 
Present work 
Present work 
Present work 
Present work 
Present Work 
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43015.8(1) 

45517.660(50) 

47221.474(20) 

48600.979(27) 

49003.629(22) 

49398.855(29) 

49797.781(22) 

50450.206(14) 

52851.33(1) 

52855.042(25) 


OSO-8 

Tenma 

Ginga 

BATSE 

BATSE 

BATSE 

BATSE 

RXTE 

RXTE 

RXTE 


Bec ker et a l.|l|1977t 
Makistiiniaet al.|fL987 


Uorbet et al. 


Kubin et al. 


Kubin et at 


Ted 


Kubin et al. 


Kubin et al. 


Uark 12UU1 


[ms 


mr 


TW7 


TW7 


TW7 


Mukheriee et al 


Baykal et al.|j20l)bf 




41449.95(7) 

Uhuru 

42628.0(1) 

Ariel-V 

43258.35(10) 

Ariel-V 

43384.889(22) 

HEAO-AI 

43563.827(22) 

HEAO-AI 

45920.2(2) 

EXOSAT 

45923.9(2) 

EXOSAT 

51016.96(3) 

RXTE/ASM 

52702.22(3) 

INTEGRAL 

52877.43(3) 

RXTE/ASM 

53779.695(50) 

INTEGRAL 

54734.18(3) 

RXTE/ASM 

54868.37(5) 

INTEGRAL 


Comii^y & Moraes||1991l 


Davison et al. 

wn 


Davison et al. 

1977 

Uommslcv & 

Aoraes 

1991 

Uommslcv & Moraes 

1991 

Uommslcv & Moraes 

1991 

Commslcy & Moraes 

1991 


Present Work 
Present Work 
Present Work 
Present Work 
Present Work 



4U 1538-522 


VelaX-l 


Epoch Time 

T 7112 

Satellite Reference 

Epoch Time 

Tnjl 

Satellite Reference 



42611.23(6) 

42727.750(24) 

42996.628(19) 

43821.34(13) 

44170.937(21) 


SAS-3 

COS-B 

COS-B 

SAS-3 

COS-B 


Rappaport et al.Hi976t 


van aer Klis & Bonnet-Bidaud ' 1984 


van Uer Ktis & B onnet- Biaau'glU984 


Kappaport et al. 


van Uer Klis Sc Bonnet-Bidaud 


1984 


44278.5466(36) 

ECHELEC 

van KerKwiiK 

et al. 

1995b| 

44305.44(7) 

Hakucho 

JNagase et al. 

1981 


45408.056(13) 

Hakucho 

JNasase et al. 

1983 


48895.2186(12) 

BATSE 

bildsten et al 

jl99 


52974.008(8) 

INTEGRAL 

Rreykenbohm et al. 

2008J 

Ted 

41446.04(7) 

Uhuru 

Forman et al. 

1973 


42449.97(2) 

Copernicus 

Charles et al. 

1978 


42620.32(3) 

Ariel-V 

Watson & Urimths 

1977 

42727.93(3) 

COS-B 

Uselman et a 

7197 

7 

43651.11(7) 

OSO-8 

Deeter et al. 

1^87 


44314.57(5) 

Hakucho 

JNagase et al. 

1983 



52176.395(16) 

52803.90(2) 

52812.88(2) 

52974.225(3) 

52983.17(3) 

53574.834(16) 

54964.333(16) 


RXTE/ASM 

INTEGRAL 

INTEGRAL 

INTEGRAL 

INTEGRAL 

RXTE/ASM 

RXTE/ASM 


Present work 
Present work 
Present work 
Present work 
Present work 
Present work 
Present work 


SMC X-1 


Epoch Time 

Tnjl 

Satellite 

Reference 




40963.99(2) 

Uhuru 

Schreier et al 

Id 19721 

42275.65(4) 

Copernicus 

luotiY & Kaplev 


975 

42836.1828(2) 

SAS-3 

Fnmini et al. 


71 

r 

42999.6567(16) 

Ariel 5 

Davison et al 

T97 

7 

43116.4448(22) 

COS-B 

lionnet-iiKlai 

S & \ 

an der Klis|jl981| 

46942.47237(15) 

Ginga 

Levine et al. 

1993 


47401.744476(7) 

Ginga 

Levine et al. 

1993 


47740.35906(3) 

Ginga 

Levine et al. 

1993 


48534.34786(35) 

ROSAT 

Woidowski et 


f 

998 


48892.4191(5) 

ROSAT 

Woi(lo\vski et al. 

mz 


49102.59109(82) 

ASCA 

Woidowski et al. 

1998 


49137.61911(5) 

ROSAT 

Woidowski et al. 

1998 


50091.170(63) 

RXTE 

Woidowski et al. 

1998 


50324.691861(8) 

RXTE 

Inametal. 2U1U 



50787.849(14) 

RXTE/ASM 

Present work 




51694.67302(1) 

RXTE 

Raichur & Paul|j2010| 

52185.052(3) 

RXTE/ASM 

Fresent Work 




52846.700(25) 

INTEGRAL 

Present Work 




52979.017(1) 

RXTE 

Raichur & PaulH2010| 

53582.268(14) 

RXTE/ASM 

Present Work 





SAX J1802.7-2017 

Epoch Time 

Tnll 

Satellite 

Reference 

52168.22(12) 

BeppoSAX 

lAugello et al.U2003| 
Hill et al. ^2005) * 

52168.26(4) 

BeppoSAX 

Ted 


52931.37(2) 

INTEGRAL 

Hill et al. 

(12005 


53260.37(7) 

INTEGRAL 

Jain et al. 

jzuuv 


53776.82(7) 

Swift 

Jain et al. 

2(109 



53863.1(14) INTEGRAL Present Work~~^ 
54503.38(7) Swift |Jain et al.||2009^ 



XTE J1855-026 

Epoch Time 

Tnjl 

Satellite 

Reference 

51495.25(2) 

RXTE 

Corbet & Mukai 2002 

Ted 

52704.038(50) 

INTEGRAL 

Present Work 

54009.972(50) 

INTEGRAL 

Present Work 
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